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Abstract 

Using the effective potential approach for composite operators, we have an- 
alytically evaluated the truly nonperturbative vacuum energy density as ob- 
tained by using a model infrared finite gluon propagator which was suggested 
by lattice simulations. The truly nonperturbative vacuum energy density is 
defined as the truly nonperturbative part of the full gluon propagator inte- 
grated over the deep infrared region (soft momentum region). With this defi- 
tion, this is a manifestly gauge invariant quantity. We have explicitly shown 
that the corresponding effective potential has always an imaginary part which 
means that the vacuum of this model is unstable. Thus this model cannot 
have a true ground state. 

PACS numbers: ll.15.Tk, 12.38.Lg 



1 



Typeset using REVT^ 



I. INTRODUCTION 



In recent papers [1,2] a general method (prescription) was formulated how to calculate 
correctly the truly nonperturbative Yang-Mills (YM) vacuum energy density (VED) of the 
ground state in QCD models by using the effective potential approach for composite oper- 
ators [3,4]. This allows one to investigate the vacuum structure by substituting some well 
justified ansatz for the full gluon propagator [| since in the absence of external sources the cor- 
responding effective potential is nothing else but the VED itself. The truly nonperturbative 
VED was defined as the truly nonperturbative part of the full gluon propagator integrated 
over the deep infrared (IR) region (soft momentum region), i.e., the method assumes that 
all kinds of perturbative contributions must be subtracted from the VED which, in general, 
is badly divergent. In order to factorize the scale dependence, the effective potential was 
introduced at a fixed scale. This makes it possible to investigate the structure of the YM 
vacuum in terms of dimensionless but physically meaningful variables and parameters. The 
nontrivial minimization procedure, which can be performed in two different ways,[| allows 
one to determine the value of the soft cutoff as a function of the corresponding scale parame- 
ter. This latter is inevitably present in any nonperturbative model (either classical, quantum 
or lattice) of the full gluon propagator. If the chosen ansatz for the full gluon propagator 
is realistic then the method will determine uniquely the truly nonperturbative VED, which 
should be always finite, negative and it should have no imaginary part (ee Ref. [1,7,8]). 
Let us emphasize that the truly nonperturbatibe VED, provided it is correctly calculated, 
should be either negative (stable vacuum configurations if it also has nontrivial minimum 
or minima) or complex (which definetely signals unstable vacuum). 

Why is it so important to calculate the truly nonperturbative VED from first principles? 
First of all, this quantity is important in its own right since it is, by definition, nothing but 
the bag constant (ie. the bag pressure) apart from the sign [9]. Through the trace anomaly 
relation [10] it assists in the correct estimation of such an important phenomenological 
nonperturbative parameter like the gluon condensate as introduced in the QCD sum rules 
approach to resonance physics [11]. Furthermore, it helps in the resolution of the U{1) 
problem [12] via the Witten-Veneziano (WV) formula for the mass of the t]' meson [13]. 
The problem here is that the topological susceptibility needed for this purpose [14,15] is 
determined by the two point correlation function from which the perturbative contributions 
should be correctly subtracted [13-18]. The same holds true for the above-mentioned bag 
constant which is a much more general quantity than the string tension since it is relevant 
for light quarks as well. Thus to correctly calculate the truly nonperturbative VED means 
to correctly understand the structure of the QCD vacuum in different models. 

In our previous publication [19] (see also the references therein) we have already inves- 
tigated the structure of the classical vacuum in the Abelian Higgs model of the dual QCD 
ground state. We have explicitly shown that the vacuum of this model without string and 

^Recent brief reviews on both the continuum and lattice gluon propagators can be found in Refs. 
[5,6]. 

^Both lead to the same numerical value of the truly nonperturbative VED provided it exists. 
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with string contributions is unstable against quantum corrections. 

Now the main purpose of this work is to investigate the YM vacuum structure by using 
an IR finite (IRF) gluon propagator as suggested by lattice simulations in Ref. [20]. Still, 
let us display a few necessary definitions, to be used later, beforehand. 



The relevant expression for the truly nonperturbative YM VED as derived (prescribed) 
in Ref. [1] is (in Euclidean metrics) 



where is the above-mentioned soft cutoff which separates the deep IR region (where the 
nonperturbative dynamics becomes dominant) from the perturbative region (see below). 
The truly nonperturbative gluon form factor is defined as follows: 



Here Aatp is the also above-mentioned scale parameter that is responsible for the nonpertur- 
bative dynamics in the model under consideration. This definition explains the difference 
between the truly nonperturbative part d^^{q^) and the full gluon propagator d{q^) which is 
nonperturbative itself. Moreover, it guarantees that the truly nonperturbative VED (2.1) is 
a manifestly gauge invariant quantity. Though the full gluon propagator is explicitly gauge 
dependent, its truly nonperturbative part is not like that since the explicit gauge dependence 
(in fact the longitudinal term which is not rcnormalizcd) vanishes along with the perturba- 
tive terms contained in (i(g^, Aatp = 0) = d^'^ (g^). (The longitudinal part, which explicitely 
depends on the gauge, is always perturbative.) It is easy to see that by "PT" we mean the 
intermediate (IM) and the ultraviolet (UV) regions (the IM region remains terra incognita 
in QCD). Fortunately the "PT" part is of no importance here. 

Thus the separation of "NP vs. PT" becomes exact because of the definition (2.2). The 
separation of "soft vs. hard" momenta also becomes exact because of the minimization 
procedure. The analysis of the truly nonperturbative VED (2.1) after the above-mentioned 
scale factorization provides, in addition, an exact criteria to distinguish between stable and 
unstable vacuum. Thus the truly nonperturbative VED, as it is given in Eq. (2.1), is 
uniquely defined. It is truly nonperturbative since it contains no perturbative information 
at all. This is rather similar to the lattice approach where, by using different "smoothing" 
techniques such as "coohng" [21], "cychng" [22], etc., it is possible to "wash out" all type 
of perturbative fluctuations and excitations of the gluon field configurations from the QCD 
vacuum in order to deal only with the true nonperturbative structure. 

The above briefly-described general method [1,2] can serve as a test of different quantum, 
classical as well as lattice models of QCD. 



Let us investigate the quantum structure of YM vacuum with the IRF behaviour of the 
full gluon propagator (in the Landau gauge) as suggested by lattice calculations in Ref. [20] 



II. THE TRULY NONPERTURBATIVE VED 




(2.1) 



rf^^(g2, A^p) = d{q\ A^p) - d{q\ A^p = 0). 



(2.2) 



III. MMS IRF GLUON PROPAGATOR 
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by Marenzoni, Martinelli and Stella (MMS). The propagator was parametrized as follows 
(from the very beginning in our notations and Euclidean metrics): 

Here M is the mass scale parameter which is responsible for the nonperturbative dynamics 
in this model, i. e., M — Anp in our notation. When the parameter M formally goes to 
zero, only the perturbative part remains. The best estimates for the parameters M and 
a are M = 160 MeV and a^^ ^ 2.0 GeV (the inverse lattice spacing). The exponent is 
7] ~ 0.53 (this mimics an anomalous dimension for the gauge field) and the renormalization 
constant is Z ^ 0.1. Let us remind the reader that within the general method [1,2], the 
nonperturbative scale parameter is considered free, i.e., as "running" (when it formally goes 
to zero then only the perturbative phase survives in the model) and its numerical value will 
be used only, provided it exists at all, at the final stage in order to numerically evaluate the 
corresponding truly nonperturbative VED. The subtraction procedure in (2.2) now looks 
like 

= M=) - M' = 0) = - ^(^.„.),i,,.^^,.(,.„.)„] . (3.2) 



A. Fixing the lattice spacing 

In order to factorize the scale dependence in (2.1), let us first choose the lattice spacing 
a as an extra scale parameter [1]. Then d^^{q^) in Eq. (3.2) becomes 

d^'^iz, h) = i, , (3.3) 

^ ' ' Zz'>[l + Zb-^z^+'i] ^ ' 

where 

z = q^a?, zq = go«^ ^ = (3-4) 

Here zq is the corresponding dimensionless soft cutoff while the parameter h has a very clear 
physical meaning, i. e., when it is zero (which means that — > 0) only the perturbative 
phase remains. 

Substituting (3.3) and (3.4) into the effective potential (2.1), one obtains 

n,{z^,h) = a^ef = --i{/i(^o,6) + h{zo,h)]. (3-5) 
where the integrals are given as follows: 

? 3 
Ii{zn,h) — dz z In [1 ; , ^ , , , ), 
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and 2-^1(15 —Zz^'^^/b) denotes the hypergeometrical function. 

From these expressions it is possible to show that the effective potential (3.5) in the 
perturbative limit {b 0) vanishes indeed, while at infinity {b — > 00) to-leading order 
becomes constant and it depends on Z , zq and i]. So at first sight, it seems the effective 
potential at a fixed scale (3.5) as function of the parameter b may have a nontrivial minimum. 

Note however that the first integral in Eqs. (3.6) always exhibits an imaginary part 
at any finite values of Zq and b. To show this it suffices to investigate the function under 
logarithm in Eqs. (3.6) 

Let us notice that R{z = 0) = —00 holds true for any fixed Z, b and 77. Since R is regular 
as a function of z in the whole open interval (0, 2:0] for any Z, b and rj^ it simply follows from 
the Boltzano-Weierstrass theorem that there exists an open interval (with z = as the left 
open end point) where R is negative, provided R becomes non-negative somewhere in the 
interval (0, Zq]. If this were not true then R must be negative in the whole interval (0, Zq]. 
Having such an interval where R < 0, and taking into consideration that the logarithm is 
a monotonous function, we certainly have an imaginary part in the effective potential (3.5) 
for any finite set of parameters Z, b, rj, zq. 

Let us derive also the formal "stationary" condition with respect to b, namely 
dflg{zQ,b)/db = 0. Then from Eqs. (3.5) and (3.6), one obtains 
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dzz^ Jl^''^]'"^,^ , =0 (3.8) 
^^{z,b)[Z<^{z,b)z'' -3] ^ ' 

with 

<f{z,b) = 1 + Zb-^z^+\ (3.9) 

From the estimates for the parameters M and a mentioned above, it follows that 6 <C 1. Still 
we use 6 < 1 in order to be sure that we do not miss the nontrivial minima. Numerically 
evaluating the "stationary" condition (3.8) (on account of the numerical values Z ^ 0.1 and 
T] f» 0.53), nevertheless we found that only the trivial solution zq = exists indeed. 

Thus the vacuum of this model is unstable and therefore it is physically not acceptable. 
The MMS IRF gluon model propagator can be related neither to quark confinement nor to 
dynamical chiral symmetry breakdown (DCSB). This statement is in complete agreement 
with the conclusion given in Ref . [23] . There it was shown that the MMS gluon propagator 
neither confines quarks nor it breaks chiral symmetry dynamically. Three special type of 
expressions for the dressed quark-gluon vertex (free from ghost contributions) were used in 
their investigation of the quark SD equation. Our result is, however, more general since we 
do not require a particular choice for the dressed quark-gluon vertex. 



^The poles at 1 + Zb ^z^~^^ = do not lie in the interval [0, zq] since zq is always positive, by 
definition. 
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B. Fixing the soft cutoff 

It is instructive to further investigate the presently discussed model by choosing the 
dimensionless variables and parameters as follows: 

-^>-i-r (3.10) 

For simplicity's sake, wc use the same notations for the dimensionless set of variables and 
parameters as in Eq. (3.4). Now when the parameter zq goes to infinity (at a fixed soft 
cutoff, see below) then only the perturbative phase survives (M — > 0). In this case, d^^{q^) 
in Eq. (3.2) becomes 

d^^(^) = (3.11) 

^ ' ai2;^[l + ai2;i+^] ^ ^ 

where 

ai = ZiMaf"^. (3.12) 
Substituting this into the Eq. (2.1) and fixing the soft cutoff itself [1], one obtains 



1 1 

^g{zo- ai) = -^ej = X ^o-2{/i(^o; Oi) + h{zQ\ Oi)}, (3.13) 

where the integrals are given as follows: 

A(-2o;ai)= dzz\n(l- . ^ 



/.(.„;a,) = j^/d.^^^. (3.14) 



A specific feature of this model is that the combination 

= Z{Maf\qlM-y = Z{qlay = v (3.15) 

is fixed when the soft momentum cutoff is fixed (like in this case, see Eq. (3.12)). Thus 
the effective potential (3.13) and corresponding integrals (3.14) become more complicated 
functions of Zq, namely 

^^9(^0, y) = = - A X Zo'{h{zo, u) + /2(^o, ^)}, (3.16) 

where integrals are given now as follows: 

? 3 
/,(.„, = y . In (1 - -.,11^ 







hyzn, v) = — — / dz = — — r2-ri 1, ; h 1; —I'Zn], (3.17 
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and 2Fi{l, —uzo) is the hypergeometrical function. 

The dependence on the parameter Zc becomes more comphcated than in Eq. (3.13) 
indeed. Nevertheless, hke in the previous case, it is possible to show again, by analysing the 
first integral in Eqs. (3.17), that the effective potential (3.16) will have an imaginary part 
(at any finite values of parameters zq and u). Consequently the vacuum arising from this 
lattice model gluon propagator is unstable indeed. 
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